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ORDER PROBLEM FOR CANONICAL SYSTEMS AND 
A CONJECTURE OF VALENT 


Abstract. We establish a sharp upper estimate for the order 
of a canonical system in terms of the Hamiltonian. This upper 
estimate becomes an equality in the case of Krein strings. As an 
application we prove a conjecture of Valent about the order of 
a certain class of Jacobi matrices with polynomial coefficients. 
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Introduction. Let L be a positive number and H be a summable function on [0, L\ 
with values in 2 x 2 matrices, such that H(x) >0 a. e. Let J — ( ^ ^ V A canonical 

system (H,L) is the matrix differential equation of the form 


(1) J— = zY\Y ; z e C, x e [0,L]. 

ax 

A solution M(x, z) of this equation satisfying M(0, z) — I is called the monodromy matrix. 
We write M(z) = M(L,z). The function H is referred to as Hamiltonian. Without loss 
of generality we assume that tr Hhr) = 1 a. e. The background on canonical systems can 
be found m [H [2]. 

Given a canonical system, the quantities 

r Dg 

hm sup- 

\z \—»oo 


\z I 


and 


hm sup 

| 2:|—>-00 


log log \ Mjj(z)\ 
log \z\ 


do not depend on i , j (see e.g. m, and are called type and order of the system, respec¬ 
tively. The type is given by the classical Krein - de Branges formula [Tj, 


( 2 ) 


type of (H,L) = f \J det H(t)dt. 
Jo 
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In particular, this formula says that if det H(T) = 0 a. e. then the matrix elements 
of M(z) have minimal type, and a fundamental question is to find or estimate the order 
of the system. This question is the order problem referred to in the title. One should 
notice here that the operators corresponding to canonical systems typically are not semi- 
bounded below hence conventional variational principles are not suitable for estimating 
their eigenvalues. 

In the present paper we establish an upper estimate for the order in terms of the 
Hamiltonian which is sharp in the power scale and gives the actual value of the order in 
all available examples where it is known. Let us formulate the result. 


Definition. A Hamitonian H is of finite rank if there exist numbers Xj, 0 = xq < x\ < 

• • ■ < xn = L, and, a finite set of vectors, ej G M 2 , of unit norm such that 

H(V) = (•, ef ) c2 e j} x G (xj, x j+1 ), j = 0,..., N - 1. 

The (elements of) sets {xj}, {ej} and the number N are called parameters and the rank 
of the Hamiltonian ,0 respectively. 


Theorem 1. Let (H,L) be a canonical system and let0<d<l. 

1. Suppose that there exists a C > 0 such that for each R large enough there exists a 
Hamiltonian H/j of a finite rank, N(R), defined on (0, L) and a set of numbers (depending 
on R) {aX iR) -\ 0 < aj < 1, for which the following conditions are satisfied (Pj = 
(■,ej)ej; x v ej are the parameters of Y\r), 

' (i) 




f* x j +1 


IIH(t) — H R (t)|| d£ < CR 


d—1 


(ll) 

(Hi) 

(iv) 


'^2a 2 j (x j+ 1 - Xj) < CR d x , 


log ( 1 + 


Pj +il 

CL j CL j _|_ 


< CR d , 


a j -1 


< CR d . 


log a 0 1 + log a Ar 1 (i?) _ 1 + log 

Then there exists a K > 0 such that 

\\M{z)\\ < e K ^ d 

for all zGC. 

2. For each p, 0 < p < 1, there exists a system (H, L) of order p which for any e > 0 
satisfies the assumption of assertion 1 with d = p + e. 


1 Notice that we do not require ej e j+i ■, hence the rank of a finite rank Hamiltonian is not defined 

uniquely. 
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Theorem 1 gives an upper bound for the order in terms of the quality of approximation 
of the Hamiltonian by piecewise constants. The choice of approximators is natural in 
the sense that finite rank Hamiltonians have order zero (the corresponding monodromy 
matrices are polynomials), see also Section [3.11 Piecewise constant (or, more generally, 
polynomial) approximations are the mainstream in studying spectral asymptotics for in¬ 
tegral and differential operators, see for instance BO El- By way of comparison, notice 
that those studies are mainly aimed at controlling the number of ’’pieces” necessary for 
approximation of a given function with a given accuracy, while in Theorem 1 the number 
N(R) does not play a direct role. In special situations, however, an optimal choice of 
approximation leads to conditions explicitly involving the number N(R) (see assumption 
(B) in the following theorem). 

An important class of canonical systems is constituted by systems with diagonal Hamil¬ 
tonians H. Such systems arise in description of mechanical strings with variable density 
sometimes called Krein strings, see [6] for details. In the context of the order prob¬ 
lem (rank H (a:) = 1) a diagonal Hamiltonian may take only two values, JV = 

and S) 2 = ll' ^ ur nex ^ resu ^ sa y s that i n this case the upper bound implied 

by Theorem 1 coincides with the actual order. The formulation is as follows. Define 
X\ = {a: G (0, L) : H(x) = T)i}, X 2 = {a: G (0,L): H(a:) = $) 2 }. Let | • | stand for the 
Lebesgue measure. 

Theorem 2. Suppose that for a. e. x G [0, L\ either H(a;) = Sji, or H(a;) = S) 2 . Then the 
order of the system (H,L) coincides with the infimum of d’s , 0 < d < 1, for which there 

exists a positive C = C(d) such that for each R large enough there exist a covering of the 


interval (0, L ) by n 

(A) 

= n(R) intervals, ojj, such that 


^2 J\uj n Xx\ |Uj nx 2 | < CR d ~ l 

(B) 



n(R) < CR d . 


We give several examples of application of Theorems 1 and 2. Namely, we estab¬ 
lish upper estimates of the order in terms of smoothness for Hamiltonians from classical 
smoothness classes (Holder, bounded variation) by applying Theorem 1, see Corollary [3J 
prove a conjecture of Valent about order of a certain class of Jacobi matrices (Corollary 
ED), and give a rather short calculation of the order for the Cantor string (see Section 1(01) . 

The first result on the order problem we are aware of is the 1939 theorem of Livsic [7] 
saying that the order of a canonical system corresponding to an indeterminate moment 
problem with moments 7 j is not less than lim sup n _ l> 00 ( 2?2 log?z)/log 72 n , the order of the 
entire function ^2 z 2j ■ A modern two-line proof of this assertion can be found in 
[ 8 ]. The next result, essentially due to Berezanskii [9], is formulated in terms of Jacobi 
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matrices. Berezanskii studied Jacobi matrices of the form 

\ 

0 ... 

P3 0 ... 

with pj growing. Although he did not explicitly address the problem of order, he made a 
crucial technical observation that allows to estimate the corresponding orthogonal polyno¬ 
mials at large j in terms of pj 1 . The explicit translation of his result to the order problem 
is given in [8]. It says essentially that if pj is a log-convex or log-concave sequence at 
large j, and qj is small relative to pj, then the order of the system coincides with the 
convergence exponent for the sequence pj. More precisely, the following assertion holds. 

Theorem 3. [8] Let |3|) be a limit-circle Jacobi matrix. If pj satisfies pj_\Pj + i < p 2 
for all j large enough, and qj/pj-i G l 1 , then the order of the system equals to inf{a > 
0: pj a G l 1 }. The same result holds if pj satisfy the inequality pj_ipj + i > p 2 instead. 

In the language of canonical systems, these results refer to a special class of Hamilto¬ 
nians defined as follows. Let bj be a bounded sequence of reals, 0 = b 0 < b\ < b 2 < ..., 
L = lim&j, and ej G M 2 , j > 1, a sequence of vectors of unit norm, ej ^ ±e 3 -_i. Let 
A j = (bj_i,bj), j > 1. Define the Hamiltonian H on (0,L) corresponding to these se¬ 
quences by 

(4) H (x) = (-,e j )e j , x G Ay 

The correspondence between Hamiltonians of this form and limit-circle Jacobi matrices is 
described in detail in [10]. Upon suitable normalization it is one-to-one, the correspond¬ 
ing selfadjoint operators are unitarily equivalent, and the Jacobi parameters qj, pj are 
expressed via ej and bj by explicit formulae. The relation of Theorem [3] and our result is 
that the relevant part of Theorem [3] (the order is not greater than the convergence expo¬ 
nent) easily follows from Theorem 1 applied to Hamiltonians of this class, see Section I3~4l 
for details. 

Apart from the mentioned general results, there are several isolated explicitly solvable 
non-trivial examples of Jacobi matrices for which the order is known and is non-zero. Two 
of them were found by Valent and his collaborators in pT] (order 1/4) and [12] (order 
1/3). In these examples, motivated by studies of the birth-death processes, q n and p 2 are 
polynomials, \q n \ ~ p n at infinity. On their basis it was conjectured in [13] that in a class 
of Jacobi matrices with polynomial p 2 and q n the order is 1/deg q n . As explained below, 
the fact that the order is not less than 1/ deg q n is almost trivial, hence the hypothesis is 
essentially that the order does not exceed 1/ degqq. We establish the latter in Corollary 
[6] applying Theorem 1. 

Another set of examples in the order problem comes from studies of uon-Weyl spectral 
asymptotics for ID differential operators. Apparently the first result in this direction is 
due to Uno and Hong |T4j who have found the order for the Cantor string. In [15], the 


(3) 


9i Pi 0 

Pi 12 P2 

0 p 2 <?3 

V 0 
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authors calculated the order (in fact, they found, in a sense, the whole leading term), for 
a class of strings with self-similar weights. The order is also known for a rather general 
class of strings related to so called d-sets [16]. Notice that nans] rely on the variational 
prinicple for the eigenvalues, hence their methods are apparently unsuitable to obtain 
results like Theorem 1 because of lack of semiboundedness. 

A general formula for the order of a string was obtained in H7j. In the situation of 
Theorem 2 it says that the order of the system (H, L) is 

{ pL /‘min {x,L—x} 

d> 0: J dM(x ) J (s(M(x + s) — M(x — s))) ^ _1 ds < oo 

Here M is a non-decreasing singular function on an interval [0,L], L + M(L) = L, such 
that X\ = {x + M(x): x G [0, L], M'{x) = 0}. This formula, to the best of our knowledge, 
has never been used to calculate the order of an actual string of the class considered in 
this paper (see also Section 16.20 . We use an argument from the proof of (J5]) in [l7j in the 
derivation of Theorem 2, see Lemma [7] 

The structure of the paper is as follows. In Section |T] we reproduce a proof of the 
inequality ’Ths of (12]) ” < ”rhs of (J2J)” from [TS] with a minor simplification. The reason we 
give it here is that it provides one of the ideas used in the proof of Theorem 1. The proofs 
of Theorems 1 and 2 occupy sections named accordingly. In the Comments sections we 
discuss the assumptions of these theorems and compare them with the earlier results. In 
the Applications section we establish upper bounds for the order in smooth classes and 
prove the Valent conjecture. 

Throughout the paper the norm signs refer to the operator norm for 2x2 matrices, 
L)i 2 are the matrices defined before Theorem 2. Unless specified otherwise summations 
extend to all values of the summation parameter for which the summand is defined. C 
stands for any constant whose exact value is of no interest for us. Given a Jacobi matrix 
(]3]), Pj( A) and Qj( A) stand for the solutions of the corresponding three-term requrrency 
relation subject to the initial conditions Pi — 1, Pq — 0, Qi — 0, Q 2 — 1/pi (orthogonal 
polynomials of the first and second kind, respectively). 



1. The upper estimate in the Krein-de Branges formula 
Proposition 1. Let (H,L) be a canonical system. Then 

type of (H, L) < f yj det H(t)dt. 

J 0 

Let p(x) be the exponential type of M(x, A). For each y e (0, L) the monodromy matrix 
satisfies the integral equation 


(6) 


Ml:,:. A) = Mill. A) - A / JHil i MU, A)d t. 
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A crude estimate of the Volterra iterations for this equation shows that | p[x) — p(y)\ < 
\x — y\ and thus p{x) is Lipschitz. The idea of the proof is to estimate p'(x) in terms of 
H and then integrate it to obtain the required bound. 

Proof. Let be a constant invertible matrix, to be chosen later. Equation (J6J) can then 
be rewritten as follows, 

f W(x, A) = QM{y, A) - A f (fiJH(f)ir 1 ) OM(f, A)d t. 

Jy 

This is a Volterra equation with respect to flM. Solving it by iterations we have (the 
Gronwall lemma), 

/ r x 

||flJH(f)f2 _1 II d t 


(7) \\QM(x, A) || < \\QM(y, A) || exp ( |A| 

for all y < x. It follows that p{x) satisfies 


p(x) < p(y) 


flJH(f)h2 d t. 


Taking the limit y f x we obtain that for a. e. x G [0, L\ 

p'(x) < 110JH(:r)Jl _1 1| . 

The Ills does not depend on 0, hence let us minimize the rhs in 
Lemma 2. Let A be a 2 x 2-matrix with tr A = 0, then 


( 8 ) 


inf 

Q : det 


IlMfr 1 !! = y/\ det A\. 


Proof. If det A ^ 0 the lemma is trivial - it suffices to choose f] to be the diagonalizer of A. 
If det A = 0, then without loss of generality one can take A = ^ ^ = diag (a -1 , a), 


and send a —» 0. 


Applying this lemma gives p'(x) < ^/det H(x) a. e., and the assertion follows by 
integrating this inequality. □ 

This proof is essentially the one in [TBJ Theorem X] except that de Branges uses an 
explicit reduction of the matrix JH(a;) rather than mere existence of a diagonalizer. 


2. Proof of Theorem 1 

2.1. The estimate. Let (H.L) be a canonical system. For an arbitrary finite set of 
numbers Xj, 0 < j < N, such that 0 = Xq < x\ < £2 < • • • < Xn = L, and arbitrary 
invertible matrices , 0 < j < N, an argument from the proof of Proposition 1 (consider 
dU) with y = Xj, x = Xj + i, Q = Llj) shows that 

( r x i +1 

|A| J JH(t)Qj 1 1| dt 
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for 0 < j < N. With the notation Mj = M(xj , A) we then have 

||fl i+1 M i+1 || < ||O i+ iQ- 1 || • ||n i M i+1 || < ||%+iQ7 1 || II QjMj 


exp | A | 


fXj+1 


fL,JH(£)fL dt 


Taking logarithm, su mmin g the resulting inequalities in j and choosing = I we get, 

jv-i , N—l 

M ° °- 1 " ' " r io| 


N 1 r x j +1 N 1 

(9) log||M(A)||<|A|£ / || QjJH(t)Qj 1 11 dt + E log 1 || + log ||U 

3=0 3=0 


Let Pj be an orthogonal rank 1 projection, Pj = (-,ej)ej, ej G M 2 , ||e_j|| = 1. Then the 
summand in the first sum in the rhs can be estimated as follows, 


r* x j+i 


f* x j+1 


( 10 ) 


||Sl 3 JH(t)nT>||d( < / ||!J,J(H(() - Pj) Up || dl + (xj+i - Xj) \\SljJPjSl- 

JXj 

r x i+i 

< im-ll Wnj'W / ||H(f) - Pj || d t + (x j+ i - Xj) \\QjJPjQj 1 


Let us choose the matrices Q r The choice is suggested by the proof of Lemma [2] 

flj = diag [aj 1 , aj) Uj , 

where Uj is a unitary transform reducing JPj into its Jordan form, 


where Uj 

and 

a j e 

/ cos 


V sin 

Vi)' 

1°. 

, ||% 

(11) 


2°. 

. For 


UjJPjUj 1 = 


0 0 
i or 


X r x j+i 

rhs of (HO]) < — / || (H(t) - Pj) || d t + a) (x j+1 - Xj). 

a j Jx , 


^3+l^j 1 — 


a j +1 0 

0 Clj -|_i 

«7+i 0 

0 dj -|_x 

a j a j +1 


‘W'T'ft ar 1 


i) J ( a i A ^ _ 


-1 - 


0 


0 dj +ia • 


-i +O 


0 a 3 


ir+i n 


a 3 a 3+l 
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Then (log(x + y) < | log 2 /| + log(l + x ) for x, y > 0) 

log||n ?+ iQj 1 || < | log (djCijU | + Clog fi + — El 

CL j CL j _|_ ^ 

Plugging this and dll]) in (J9]) and taking into account that ||f2o|| = % 1 1 H^tv-iII = <%-i> 
we obtain the first assertion of the theorem. 

2.2. Sharpness. Let p G (0,1), a = p _1 — 1, d — p + e. Define bj — 1 — j~ a for j > 1, 
and for x G [0,1] let 

-fin x ^ Uj&j-i, b2j] 

S) 2 , x ^ Uj[^2j—l, b 2 j\. 

The required assertion will be proved if we show that (a) H satisfies the conditions 
(i)-(iv) of the theorem for all e > 0, (b) the order of the system (H, 1) is not less than p. 
Given an R > 0. define 

H(x), x G [0,6jv-i] 

•Qi, x G [6jv_i, 1], 


H(a;) = 


Hij(x) = 


with N = N(R) to be chosen later. Let a^-i — 1, aj — R 2 for j < N — 2. We then 
have, 

lh80f(l) 5 (n - ir = 0 PT ■ 

N—2 

lhs of (h) < 4r 1_ “ + N ~ a = 0 ( Rd ^) + 0 ( N ~ a ), 

0 

lhs of (iii) < C(N - 1) log R + 0(log R) = O (N log R ), 
lhs of (iv) = 0(log R). 

Let N ~ R p as R —> 00 . Then assumptions (i)-(iv) are satisfied. 

Let us now establish that the order of the system is not less than p. To this end, we 
use the following identity. Let 0(x, A) be the first column of M(x, A). Differentiating 
(0, J@) c2 with respect to the equation dH) , we find 

(12) s(M n (l,A)M 21 (l,A)) = 3A f 1 <H(*)0(t,A),0(t,A)> c2 di. 


In the situation under consideration the rhs is 

c> A Sr (h h \ M u (bj,\), j even 

M 21 (bj, A), j odd 

For 3 A > 0 one can estimate this quantity from below. Let S 3 = bj — bj_ 1 . We have 
(13) rhs of (fill) > SA S 2 j \Mn(b 2 j, A)| 2 = 3A 5 2 j \Mn{b 2 j-i, A)| 2 . 
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In the last equality we took into account that Mu ( 627 , A) = Mu ( 627-1 j A) in the situation 
under consideration as M[fx, A) = 0 when x G ( 627 - 1 , 6 2 j). 

To estimate the rhs of (1T3T) from below we use the following corollary of the fact that 
all the zeroes of matrix elements of M(x, A) are real. 

Remark 1. Let ( G,L) be a canonical system such that (0, L) is a union of disjoint 
intervals, Ij, accumulating only at L, and G(x ) is a constant rank 1 operator on each 
Ij. Then M(x, •) is a polynomial for all x G (0, L). For any m,l, 1 < m,l < 2, define 
k{x) to be the degree of the polynomial M m i(x, ■), c(x) be its leading coefficient, 

M m i(x , A) = c(x)X k ^ + (a polynomial of degree < k(x) — 1). 

Then \M m i(x,ir)\ > |c(x)|r fc ^) when r > 1. 


Applied in the situation under consideration to the left upper entry of the monodromy 
matrix at x — 627 - 1 , this gives |Mu( 627 - 1 , ir)\ > \cj\r kj for r > 1 , c 3 and kj being the 
leading coefficient and the degree of the polynomial Mu( 627 - 1 , •), resp., j > 1 . 

Let us calculate kj and Cj. Define Mj( A) to be the value of the matrix solution of (JT|) 
with the Cauchy data Y(bj) = I at x = 6 J+i , then 

(14) M(6 2i _i,A) = M 2i _ 2 (A)M 2 7_ 3 (A)---M 2 (A)Mi(A) 


by the multiplicative property of the monodromy matrices. A straightforward calculation 
gives 


M,(\) 


I + A6j_|_i 


0 0 

-1 0 

0 1 

0 0 


j odd 
j even. 


The leading term in the matrix polynomial M( 627 - 1 , A) comes from choosing the terms of 
the first order in A in each multiple in (1141) . It has the form 


(-1) J+1 <%—1627—2---6 2 



Thus, kj = 2 j — 2, \cj\ = \^=2 an d one can continue the inequality in (fT3]h 


^27-1 


rhs of (H3D > ^3 n « IAI 


2(27-2) 


, n =2 


On the other hand, if p is the order of the system, then for any e > 0 the lhs in (fT21) 
is not greater than exp ( C £ r p+e ). By a standard relation between Taylor coefficients and 
exponential order (s. f. [19]) it follows that for large j 


J 27 


f2j-l 

lit 

, 71=2 


4 j 

<l f^ 
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Notice that the sequence 8j = j~ a — (j + l) -a is monotone decreasing in j, hence the last 
inequality implies that 82 j = O (jD 1,/(p+e ' 1 ). Comparing this with 8 j x j~ 1 ~ a = j~ 1//p , we 
hnd p + £ > p for all e > 0, that is, p > p. The proof of Theorem 1 is thus completed. 


3. Comments on Theorem 1 


3.1. Choice of approximants. Let us consider a canonical system having Hamiltonian 


of the form H = {-,e(x))e(x) with e{x) = 
subject to the condition u'v — 
Hamiltonian is 1/2, for if Y = 


u{x) 

v{x) 


, u,v being smooth functions on (0, L) 


subject to the condition u'v — v'u = —1. The order of the canonical system with this 

'Y 

+ ' is a solution of the system, then a straightforward 


Y_ 


calculation J20] shows that y = Y + u+Y_v satisfies the Schrodinger eqnaiton —y"+qy = A y 
with the potential q = u"/u. This suggests that smooth functions cannot be used as 
approximants, at least in the whole range ( 0 , 1 ) of orders. 


3.2. Formulation. The Krein - de Branges formula implies that if assumptions of The¬ 
orem 1 are satisfied for some d < 1 then rank H (x) = 1 a. e. This fact can eas¬ 
ily be seen directly. Indeed, suppose that (H,L) is a canonical system such that for 
any R large enough the conditions (i)-(iv) are satisfied for some d £ (0,1). Define 
S e = {t: ||H(t)/|| > e||/|| for all / £ C 2 }, e > 0. Arguing by contradiction, let e > 0 be 
such that |5 e | > 0. Applying the Schwarz inequality and using conditions (i) and (ii) we 

hnd that Y \J{ x j +1 ~ x j) f ^ +1 ||H(t) — H^(t)|| dt < CR d ~ l . The quantity \[xj,Xj + i] D S e \ 

estimates from below both factors in the summand, hence the sum is not less than 
\\ x ji x j+ 1 ] Ll S e | = |5 e | > 0. Taking the limit R —> 00 we obtain a contradiction. 

3.3. Sharpness. In the example establishing part 2 of Theorem 1 the conditions (i)-(iv) 
are satisfied with d = p if we insert the log A factor in the rhs. It is not known to the 
author if one can get rid of the logarithmic factor, that is, part 2 holds with e = 0 . 

3.4. Comparison with Theorem [3l The relevant part of Theorem 3 is the assertion 
that under the stated assumptions if {pj 1 } € /“ then the order < a. Let us first translate 
the setup of Theorem 3 to the language of canonical systems. As mentioned in the intro¬ 
duction, given a Jacobi matrix (l3lh explicit formulae that express the Jacobi parameters 
qj, pj via the corresponding Hamiltonian are known, see [10]. Let us reproduce them in a 
convenient form. In the following theorem P n and Q n are the orthogonal polynomials of 
first and second kind defined by the matrix (J3J) , respectively, and e 3 and bj are parameters 
of the corresponding Hamiltonian of the type described after Theorem 3. We write ipj for 
the argument of the vector Cj £ M 2 , A j = (bj-i,bj), 8 j = |Aj|. 

Theorem 4. [[TDj The correspondence between canonical systems and limit-circle Jacobi 
matrices can be chosen so that 
( 1 ) 8 n = P n (0) 2 + Q n (0) 2 , 
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(ii) 


Pj 


sin 


1 

fj+ 1) I y/fij+ifij 


3> I- 


Now let qj = 0 for simplicity. Then it can be shown [10] that Cj _!_ Cj-i, e,\ — ^j, 
hence pf 1 = y/SjSj + i. Solving for Sj we have 




Pj-lPj-3 


Pj Pj ~ 2 

Let pj satisfy the assumption of Theorem [3] Then Pj-i/Pj at large j is a monotone 
sequence having a limit < 1, and if the limit is 1, then it is increasing. This implies that 
5j = O (pjl i) and therefore {(/,} G l d for any d greater than the convergence exponent 
of the sequence pj. Notice that the mondromy matrix corresponding to the interval A j 
is I + O (|A|<Sj). This, the multiplicative property of monodromy matrices and {</,} G l d 
imply by elementary inequalities that the monodromy matrix of the system in question 

is O (e c 'l A l J, which gives the assertion under consideration. 


It is just as easy to derive this assertion from Theorem 1. Let d be such that {p ? x } G l d 
and let 01 = {j: Sj > R -1 }. Define Hr(x) = H(a;) whenever x G A*, k G 01. On 
the complement of we define to be an arbitrary constant rank 1 orthogonal 

projection. With this definition, is a finite rank Hamiltonian with parameters to be 
denoted Xj , /j. Let aj = whenever j is such that [xj,Xj + 1 ] coincides with one of 

the intervals A^, k G 01, aj — 1 otherwise. Then conditions (i) and (ii) in Theorem 1 are 
satisfied because {<5j} G l d and so Ylsj<R-i^j = ® conditions (iii) and (iv) are 

satisfied because the number of f s for which Sj > R 1 is O ( R d ), again by {<5j} G l d , and 
therefore the rank of H/j is O ( R d ) as well. Applying Theorem 1 we conclude that under 
the assumptions of Theorem [3] with qj = 0 the order is not greater than d. Thus, our 
result generalizes the upper estimate in Theorem [3] 

The case of qj subject to the smallness condition of Theorem [3] can be obtained from 
this by standard methods of abstract perturbation theory. 


4. Applications 

4.1. Smooth classes. A corollary of Theorem 1 is obtained when the conditions (i)-(iv) 
are satisfied with aj independent of j. In this case condition (ii) reduces to a? = O (i? d_1 ) 
and the lhs in the other three conditions is monotone decreasing in aj, hence without 
loss of generality one can assume that all four conditions are satisfied with aj = R d ~ 1 . 
Introduce the following notation. Given a finite rank Hamiltonian, G, with parameters 
{ej, let Var G : = \\P 3 - P j+ i||, Pj = (•, e,)e r 

Corollary 3. Assume 1/2 < d < 1. Let for any e > 0 a finite rank Hamiltonian H e 
defined on (0, L) exist such that 
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(a) 

IIH — H IIl i (o,l) — e -> 

(b) 

Var H e < Ce 53=5. 

Then the order of the system (H, L) is not greater than d. 

Proof. System (H, L) obeys the condition of Theorem 1 with H# = H 6 ^, e(R) = P 2 ( d_1 ), 
aj = Condition (a) implies (i), (b) implies (iii) via an elementary inequality, (ii) 

and (iv) are immediate. □ 

This corollary allows to give an upper estimate for order of Hamiltonians in classical 
smoothness classes. We give two examples. 

Corollary 4. Let (H,L) be a canonical system with rankH(x) = 1 a. e. 

1. If H G C“[0, L\, 0 < a < 1, then the order of the system is not greater than 1 — a/2. 

2. If H has bounded variation then the order is not greater than 1/2. 

Proof. The first assertion follows from choosing Xj = Lj/N, ej G Ran H(xj), for the 
parameters of the approximating Hamiltonian H e and adjusting N. The second assertion 
is trivial. □ 

Notice that while the estimate of Theorem 1 is sharp, the Hamiltonian in the corre¬ 
sponding example is discontinuous. It is an open question whether the first assertion of 
Corollary [I] is sharp. The second assertion of the corollary admits an ’’elementary” proof 
based on a trick from |2Tj, Theorem 3.6]. 

4.2. Berg—Valent matrix. In this subsection we consider the order 1/4 Jacobi matrix of 
Berg-Valent HU as a warmup for the proof of the Valent conjecture. We do not need the 
explicit formulae for the parameters q n and p n . The necessary information about them 
from HU is as follows 
1 °. p n ~ n 4 as n —* oo. 

2°. The values of the corresponding orthogonal polynomials, P n (z) and Q n (z), at z = 0 
have asymptotics P n (0) ~ epu -1 , Q n (0) ~ C2UG 1 with c \ y 2 7^ 0 [II] (2.33), (3.2)]. 

By Theorem H] we find that 5j ~ Cj~ 2 , sin {tpj — tpj+f) = 0(j~ 2 ). Let H# = H on 
(0, b N _i ) with TV ^ 1/4 Th d Ti define arbitrarily on (6/v—1, L) so that 

H/j becomes a finite rank Hamiltonian on (0, L). Let etjv-i = 1- Then 

lhs of (i) < 2 (L - 6 jv-i) = 2 = 0 (P d_1 ) , 

j>N 

N ~ 2 a 2 

lhsof(ii) = + ° (^ _1 ) , 

1=0 3 

lhs of (iii) < ^ log f 1 + 0 - 

- ^ V J 2 a jaj+l 


+ 0 ( 1 ). 
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Define a 2 = R d 1 for j < R d , a 2 = R 2d 1 for R d < j < N — 2. Then Y2o 2 a 2 j 2 = 
0 the lhs in (iv) is 0(logi2), and 

^ log (l + ^) < CR d log R + R l ~ 2d Y,~ 2 =° ( Rd lo § R ) + 0 ( Rl ~ 3d ) = 

V J a j a j+iJ j>Rd J 

O {R d log /?) 

because d > 1/4. Applying Theorem 1 we conclude that the order is not greater than 
1/4, the actual order of the system found in [TTj . 

4.3. Valent’s conjecture. The following assertion generalizes the consideration of the 
previous example. 

Proposition 5. Assume that a Jacobi matrix (Q|) is such that P 2 { 0) + <3 2 (0) ~ Cn A ~ D , 
p n ~ n D as n —> oo with numbers A, D satisfying 1 < A < D — 1. T/ien t/ie order zs not 
greater than 1/D. 

Let A n ,/j, n , n > 0, be sequences of reals, A„ > 0 for n > 0, p n > 0 for n > 1, p 0 = 0. 
Define 

(15) Qn +1 A n T Pm Pn+1 \/A n /X n _|-i. 

The Jacobi matrix with parameters qj, pj is said to be corresponding to birth-death 
processes with rates X n and p n m- 

Corollary 6. The order of the Jacobi matrix corresponding to birth-death processes with 
polynomial rates X n = (n + Bf) ■ ■ ■ (n + Bf), p n = (■n + Af) ■ ■ ■ (n + Af) subject to the 
condition 1 < ^f{Bj — Aj) < £ — l, is 1/ 1. 

Let us establish the proposition first. 

Proof. By Theorem [4] the assumption of the proposition implies that 

~ Cj A ~ D , sin fa - <p j+1 ) = O (j~ A ) . 

Fix a d > D -1 small enough and define as in Section 14.21 The value of N is to be 
chosen so that Y2j>N Sj x jyA-D+i _ q ( R d- thus j e |. y ^ _Ra-d+i . Define 

r R d -\ j < R d , 

a 2 = < R d—l+d(D—A—l ), R d < j < N __ 1 

(l, j = N- 1 

Notice that R d A — 1, so the corresponding range of j’s is non-empty. With this choice, 
(i) of Theorem 1 is satisfied by the choice of A, (ii) is satisfied because ^/?rf<j<v_i = 
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O (R d ( A D+1 )) and aj for this range are chosen precisely to make the corresponding term 
O (i? d_1 ), and the ills in (iv) is O (log R). The lhs in (iii) is estimated above by 


lo s ( x + - 




< CR d logi? + 


ftd-l+d(D-A-l) x_, j 


£ is = O (R d+ ’) + 

j>R d ' 

O (i? 1_dD+<i ) . 


The rhs is O (/? d+£ ) for d > D 1 and any £ > 0, and the assertion of the proposition 
follows by Theorem 1. □ 


Proof of Corollary The fact that the order does not exceed 1 ji follows from Propo¬ 
sition [5] by inspection of (fT5]l and explicit formulae [TTj expressing Pj(0 ), Qj(0 ) via Aj’s 
and /ij’s. On the other hand, an application of j8, Propositions 7.1] shows that the order 
is not less than 1 ji. For completeness, we provide a proof of the latter fact. First, for 
any p greater than the order of the system there exists a K > 0 such that for all £ large 
enough 

£ WMI 2 < e*W P . 

This is an easy corollary of the Kristoffel-Darboux formula. Since all zeroes of Pf s are 
real, \Pj(ir)\ > TTjC for r > 0, 7Tj = 1/ (pi • • • Pj) being the leading coefficient of the 
polynomial Pj, therefore 

( x -i 

i) 

c) 

Under the assumptions of the corollary, pj ~ j l , which implies p > \jl. □ 

The assertion of Corollary [6] was conjectured in [13] on the basis of two explicitly 
solvable examples, the one dealt with in the previous subsection and another one [12] 
with t = 3. 


5. Proof of Theorem [2 

The structure of the proof is as follows. First we are going to show that the order of the 
system is not greater than the inhmum. This will be done by an application of Theorem 
1 to a natural approximation H^. Then we will show that the order is not less than the 
inhmum by an appropriate choice of the covering. 

The order < the infimum. Let d be such that for some C > 0 for each R large enough 
there exists a covering of the interval (0, L) by n(R) < CR d / log R intervals, to be denoted 
ujj, such that (A) is satisfied. The stated inequality will be established if we show that the 
order is not greater than d. Without loss of generality on can assume that the intervals 
ojj are mutually disjoint. Define 

(<di, x E Uj, \uj fl Xl| > , 

1 5)2 otherwise. 
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With this choice of H 


R 


|| H(t) — Hft(t)|| dt < 2 minjlcUj D Xx\, \ujj fl X 2 |}, 


hence the condition (i) of Theorem 1 takes the form 

1 


^ — min{|a;j fl Xi\, \ujj fl X 2 \} < CR 


d -1 


Notice that min{|o;jnX l |, lo^flA^I} x \ \ujj nA 2 |/|u;j|, hence the latter condition 

is equivalent to 

1 |tm n X 1 1 I Uj fU 2 | < 


(16) 


E 




Conditions (iii) and (iv) of Theorem 1 in the situation under consideration are satisfied if 
(17) ^logfl+oj 1 ) <CR d , 

and condition (ii) has the form 

(is) J2 a2 j\^\ - 

Let a,j = 1 for \ujjl < 2/R. We write 91= {j : \ujj\ < 2 /R}. The parts of sums in 
(USD, (fT71) and (1T8T) over j G 91 are then estimated above by n(R)R ~ l , n(R)R~ 1 and n(R), 
respectively, hence they are O by condition (B). For \ujj\ > 2/R we optimize the 

choice of a 3 over the summands in (1T6|) . (1171) and (fTST) by taking 


dj = max 


j l \f\wj n Ai 

| ujj n x 2 \ | 

1 R\ujY 

\Uj 

i i 


With this choice the sums over j ^ 91 in (USD? (HID and (fl7|) are estimated above by 


'y ] J|Uj n X \| \ojj n x 2 \, 


R-'niR) + E\A ojj n x 111 ujj n x 2 1, 

and 

^log^Wjl) < n(R) log R, 

respectively. Plugging here (A) and (B) of Theorem 2 we obtain that all the assumptions 
of Theorem 1 are satisfied. 

The infimum < the order. Given ar>0,i6 (0, L), let s(t,x) G [0, a;] be the solution 
of the equation r 2 |(s,a:) fl Aj | |(s,ic) fl X 2 | = 1. This solution is unique whenever exists. 
Without loss of generality one can assume that, say, for some a > 0 the function H(x) = fp 
for x G (0,a/2), H(a;) = S) 2 for x £ (a/2, a) (attaching such two intervals at the left end 
does not change the order). Then s(r, x) is defined for r large enough for all x > a. 
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Lemma 7. |jI7[ lemmas 1-3] The order of the system (H,L) equals to 

( L Xl(x) 


(19) 


lim sup 

r—>•+oo 


|(s(t,x),x) fl X 2 
l°g T 


-dx 


X 2 being the indicator function of the set X 2 . 

This assertion provides a crucial step in the proof of the Rats’ formula for the order. It 
is formulated in in terms of the corresponding strings. For completeness we reproduce 
here the proof of the part of the lemma that we are going to use - the order is not less 
than the quantity (fT9lh translated to the language of canonical systems. 

Proof. By definition the order of a system is the order of any of the matrix elements of the 
monodromy matrix. Let us show that the order of the matrix element M n (z) is estimated 
from below by the rhs in (fTUjh The order of M u coincides with 

log log | Mn (ir) | 


lim sup 

T—>-+00 


log T 


because Mu is a real entire function having all its zeroes real. 

Let Xi — 1 — Xii p( s ,t) ~ |(s,t) Fl X 2 \. On rewriting the first column of ([I]) as an 
integral equation we obtain, 

rx nx 

M u (z,x) = l + z X 2 (t)M 21 {z,t)dt = 1 - z 2 Xi(t)p(t,x)M u (z,t)dt. 

Jo Jo 

When z = ir, t > 0, this becomes (£ r (x) : = Mn{ir,x)), 


£ r (x) = 1 + r- 


Xi(t)p(t,x)£ T (t)dt. 


It shows that is a positive and monotone nondecreasing function. Let us estimate 
£'(x)/£ T (x) from below. For a. e. x G X 2 and all s < x, we have 

/ S PX 

Xl(t)Crdt~T 2 J p(t,x)xi(t)£r(t)dt 

/ .S PX 

Xlfrdt - T 2 p(s , x) J Xltrdt = &(x) ~ p(s, x)£(x), 

and thus 

( px ps\ px 

J + J ) Xlfrd t > T 2 J Xlirdt > T 2 |(s, X) fl X ± \ f T (s) 

> T 2 |(s,x) nXi| (^t(s) -p(s,x)£(x)) . 

Picking s = s(r, x) we obtain that for a. e. x G X 2 fl (a, L ) 


£-0*0 

£0*0 “ 2 p(s(t,x),xY 


> 
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which implies the required assertion upon integration in x over [a, L\. 


□ 


( 20 ) 


The proof of Theorem 2 will be completed if we show that for any d such that 

f L X2(x) 


-dx = O (r d ) , r —y +oo, 


Ja P(S(T,X),X) 

the interval (0, L) can be covered by O (f? d log R) intervals, Uj, so that 


( 21 ) v/l Uj n Xi| I Uj n x 2 \ = O (R ^ 1 ). 

For each R large enough define a monotone decreasing sequence {x 2 }, j > 1, as follows, 
x\ = L, Xj+ 1 = s(R,Xj), if j > 1 and Xj > a; if Xj- 1 > a, Xj < a then x 3+ \ = 0 and 
the sequence terminates. Observe that the sequence x 3 is finite. This follows from the 
definition of the function s(R,x ), for 

1 = R 2 \ (s(R, x), x) fl Xl| |(s(i2, x),x) D X 2 \ < R 2 |x — s(R, x) | 2 , 

which means that x 3 — x 3+ \ > R so the sequence has O(R) members. Dehne ujj = 
[x ]+ \, x 3 ]. By construction [0, L\ = U jUJj. We claim that this is the required covering. 
First, we have to show that N(R), the number of intervals in the covering, is O (yR d logi?). 
To this end, notice that p(s(r,x),x) < p(xj + 2 ,Xj) for x G [x J+] , x 3 \ hence 


N(R)-1 

lhs in d2HD > Y' -—-, 

j = i s j + s j +1 


3 


= p{x j+X ,Xj). 


Let Q = {j: < 2}, g = {j: > 2}; n s , n g be the respective numbers of elements. 

When j G g the summand in the displayed sum is bounded below, hence, n g is 0(R d ) 
by (gUP - To estimate h g notice that s 3 > 1/ ( LR 2 ), for 1 = R 2 \uj 0^1 s 3 < R 2 Ls 3 . It 
follows that if k is the length of a discrete interval of the set g and m is the right end of 
it then x m — x m+ i > 2 k ~ l / (LR 2 ). On the other hand, x m — x m +i < L trivially, hence 
k < C + 2 log 2 R, and therefore h g = O ( R d log R). Thus, N(R) — n g + n g — O ( R d log R) 
as required. To complete the proof, notice that by the very definition of s(R, x ), the 
summand in (l2Tjl equals to hence (|2TT) reduces to N(R) = O (i? d ) and thus holds 
trivially. □ 


6. Comments on Theorem 2 and applications 


6.1. The Cantor string. Let £: [0,1] —» [0,1] be the standard Cantor function, (■ Uj,Vj) 
be its constancy intervals, T(x) = x + £(x), L = 2. Dehne the Hamiltonian H on [0, 2] to 
be 

3i> jT^u^Vj]) 

$) 2 otherwise. 


Hfx) = 


The canonical system (H,2) is called the Cantor string. We are going to show that the 
assumptions of Theorem [2] hold for d = dc = 2/ log 2 6. Let t 3 be the union of 2 J ^ 1 intervals 
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thrown away on the j-tfi step of construction of the Cantor set. Define Mr = T (ULi r k ) , 

j ~ d log 2 R- The set Mr is a union of 0( 2*) non-intersecting intervals. Consider the 
covering of [0, 2] by the intervals of the set Mr and their contiguency intervals, the latter 
to be denoted c jj. By construction, the overall number of intervals in this covering is 

0 ( 2 *') = O ( R d ). 

The terms corresponding to the intervals of the set Mr in the sum in condition (A) 
obviously vanish, hence the sum reduces to 

\J|Wj Cl Ad| |ujj Cl X 2 | < \ujj PI Aj| = \/|((0, 2) \ Mr) D Ad| 

Since T is linear on intervals of sets t*,, T~ 1 A" 1 = t v and therefore 

I((0,2) \ Mr) n Adi = \n\ = Qy ~ = R 2 (d- 1 ) 

k>j 

It follows that the assumption (A) is satisifed. Applying Theorem [2] we conclude that 
the order of this system is not greater than dc ■ Let us show that the order is not less 
than dc- Assume that the condition of Theorem 2 is satisfied for some d and let us use 
D = = {Aj} for the covering whose existence is required by the condition. Let 

j = d log 2 R and Mr be defined as above. Then without loss of generality one can assume 
that the intervals A j are mutually disjoint and the intervals of Mr are among them. 
Indeed, adding the intervals of Mr to D and removing the intersections of intervals of T) 
with the intervals of Mr does not increase the lhs in (A) and increases the constant in the 
rhs of (B) by at most 1. Given an u, an interval of contiguity of Mr, let n w ,R be the number 
of intervals of D belonging to interval oj, then n Wj R < CR d = C 2 J . Since the number of 
intervals of contiguity is 2*, it follows that the number of fc’s for which tt^r < 2 C is greater 
than 2* _1 . Let w be a contiguity interval for which nk,R < 2C, jo = j + log 2 (2C) + 4. 
By a Dirchlet box argument then there exists an interval A 6 3 which contains two 
nearby intervals of T (rj 0 ) , for the number of intervals of T (r JO ) contained in cu is 2* 0- *. 
By construction the measure |A D Aj| > 2 • 3 - * 0 = C3 ~ j , and |A fi X 2 \ > 2~ jo = C2~\ 
so yj\ A fl Xi| | A D Ad| > C2~^ 2 3~^ 2 . Since the number of such intervals A is estimated 
below by C2* we find that the lhs in (A) is estimated below by C (2/3) i/2 = CR d ^- l °^ 3 ^ 2 . 
This being O (A d_1 ) means that d > dc- Thus we have derived that the order equals to 
dc, the result obtained in [14] or [T5J by other means. 

6.2. Kats formula. As mentioned in the introduction a direct comparison of Theorem 2 
and Kats formula (|SJ) is not possible for lack of examples using the latter in the situation 
of the former. Notice however that, properly understood, (JHj) holds for a class of problems 
with L — oo (singular strings) and examples are known [22] in this class where the order 
is calculated by application of flH]). 

2 The factor 2 in the numerator in the expression for dc is due to the spectral parameter in our definition 

of the string being the square root of a ” natural” parameter used in m- 
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